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We study the effects of adding loops to a critical percolation cluster on the dif- 
fusional, and equivalently, (scalar) elastic properties of the fractal network. From 
the numerical calculations of the eigenspectrum of the transition probability matrix, 
Q\ . we find that the spectral dimension d s and the walk dimension d w change suddenly 

as soon as the floppy ends of a critical percolation cluster are connected together to 
form relatively large loops, and that the additional inclusion of successively smaller 
loops only change these exponents little if at all. This suggests that there is a new 
universality class associated with the loop-enhanced percolation problem. 
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I. Introduction 

One of the few experimentally accessible measures of the degree to which a struc- 
ture is fractal 0, § is the spectral dimension d s |J. The spectral dimension was 
originally defined 0] from the low energy behavior of the vibrational density of states 
n(E) of a fractal, elastic network: 

n(E) ~ E ds l 2 ~ l , (1) 

where E is the energy of a mode. The energy is proportional to u> 2 where uj is the 
angular frequency of the vibrational mode, and writing the above asymptotic relation 
in terms of the density of states per unit interval in u, we have the relation 

PH ~ w*" 1 , (2) 

in which d s has replaced the Euclidean dimension d in a more familiar equation for 
the phonon density of states. 

The exponent d s can be measured by various optical and neutron scattering ex- 
periments H , and for some silica aerogels in particular, a value in the neighborhood 
1.3 ± 0.01 was found. Now, such a fractal is sometimes compared to the critical 
percolation cluster || since the latter is the simplest and best studied model of an 
equilibrium random fractal. If this comparison is made for the silica aerogel, the 
experimental value of d s happens to be close to the d s of the scalar elastic problem J3J 
on the critical percolation cluster in d = 3, and much larger than the corresponding 
vector elastic result of about d s = 0.9 [j]. (Vector elasticity is the problem where the 
displacement of each node in the network is a vector rather than a scalar.) Courtens 
H argued that this is not an evidence for scalar elasticity but rather that an aerogel 
has very few floppy endings unlike a percolation cluster. That is, most floppy endings 
which may form during the growth of an aerogel are connected together to form loops 
during the aging and (supercritical) drying processes. Thus the final structure is much 
more rigid with higher value of d s although the fractal dimension df is essentially the 
same. 

In this paper, we test this suggestion, but for the simpler, scalar elasticity problem, 
on two- and three-dimensional critical percolation clusters. In this case, each bond 
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connecting two (occupied) nearest neighbor sites of the cluster acts like an ideal spring 
of uniform spring constant in the sense that it exerts a Hooke's law force based on 
the difference in the (scalar) displacements of the two sites. The effects of loops are 
studied by adding controlled amount (and size) of loops and calculating the dynamic 
exponents d s and d w @. As far as we know, this is the first quantitative study of the 
effects of adding loops on the structural properties of a stochastic fractal. 

We study this problem as a diffusion problem using the vibration-diffusion anal- 
ogy H for scalar elasticity. According to this correspondence, e.g., the vibrational 
spectrum of the network can be calculated as the spectrum of the transition probabil- 
ity matrix W M, HL where the elements Wij of W is simply the hopping probability 
per step from node j to node i in the discrete network. Thus this matrix W contains 
the information on both the network topology and the dynamics of the model of dif- 
fusion. We choose for our calculations the so-called blind ant random walk || as a 
model of diffusion although the specific choice of the type of random walk is irrelevant 
for our results. For a blind ant, the random walker attempts to move without the 
knowledge of which neighbors are part of the cluster (thus available) and which ones 
are not, and thus when it happens to choose an unavailable neighbor, it cannot hop 
and must wait at its current position for the next time step. 

From the point of view of the diffusion problem, the interesting quantities include 
the probability P(t) for a random walker to return to its starting point after t steps 
and the root-mean-square distance R(t) traveled by the random walk in time t. These 
quantities are expected to have the power-laws for asymptotically long times t: 

P(t) ~ r ds/2 (3) 
R(t) ~ t 1/dw . (4) 

These relations may be considered to define the exponents d s and d w . 

We calculate these exponents by approximately diagonalizing the matrix W using 
the method of Ref. ||, §] . Once the diagonalization is done, we compute two quan- 
tities, the density of eigenvalues n(A) (where A denotes the eigenvalues of W) and a 
certain function 7r(A) (which is the product of n(\) and some coefficient determined 
when the stationary initial state distribution is expanded in terms of the eigenvectors 
of W PT0|1 ) . These functions are expected to behave, asymptotically near A = 1 [|J, 
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as 



n(A) ~ (lnAI^/ 2 - 1 (5) 
tt(A) ~ llnAI 1 " 2 ^. (6) 



II. Loop Addition 



In this section, we characterize the clusters for which the calculations described 
in Section I are performed. The procedure is as follows: we first generate random 
realizations of occupied sites with the known critical probability p c (independently 
of each other) on either a square or simple-cubic grid of a predetermined size (edge 
length L). The percolation problem is then constructed by imposing the nearest- 
neighbor connectivity among the occupied sites using periodic boundary condition in 
all directions. We then choose only those realizations whose maximal cluster spans 
the grid in all directions as well as wraps in all directions, i.e., that those on which 
there is a path that winds around the grid (with the periodic boundaries) in each 
coordinate direction. The maximal, wrapping cluster generated in this way is an 
allowable starting configuration. 

In order to add loops in a controlled manner, we first mark all perimeter sites 
which are neighbors to two or more occupied sites (i.e., those empty sites which 
would connect two or more occupied sites if they were occupied). For each of these 
multiple-perimeter sites, we must decide the size of the loop it closes. While it is very 
easy to determine the size of the smallest loop, that is not a very interesting quantity, 
since it may also close a much larger loop and such a larger loop may dominate the 
rigidity of the structure. On the other hand, the size of the largest loop itself is also 
not an interesting quantity since in principle there are many loops that are embedded 
in the interior of the cluster with little influence on the cluster's rigidity. 

Thus we use the following procedure to determine the effective size of the loop a 
multiple-perimeter site closes: For each such site, we consider every pair of occupied 
neighbors in turn. For each pair, we calculate the minimum loop size by a burning 
method [iT|, spreading a fire from each end and determining the first contact of the 



two fires. Then we compute the maximum among these minimum loops sizes attached 
to the perimeter site in question. This is the appropriate loop size P, the size of the 



floppiest loop associated with the particular multiple-perimeter site. We then add all 
multiple-perimeter sites with P greater than or equal to a certain predetermined value 
P Q to control how floppy a loop must be for it to be closed. For example, P D = oo 
will close no loop, and P = 4 will close all loops for both the square and simple cubic 
lattices. 

In this paper, we compare the results for P a = oo, 12, 8, and 6. A typical starting 
cluster and the corresponding loop-added clusters on the square lattice are shown in 
Fig.l. Clearly, the fuzzy and floppy endings present in the original cluster (a) become 
thicker and better connected successively in (b), (c), and (d) as more loops are added, 
giving the impression of image sharpening. The final structure in (d) with P a = 6 
appears to be a much more solid object than the original; yet, the fractal dimension 
of all these clusters are virtually the same, being about 1.9. 

The above procedure is followed after the starting configuration is fixed. The 
applicable multiple-perimeter sites are ordered in a particular way (in a typewriter 
fashion due to the construction of the cluster), and the addition procedure applied 
only once in the fixed order, with no recursion. That is, once a site is determined 
to be either added or not, we do not later consider it again for addition and also we 
do not consider perimeter sites newly created by this process itself. This removes 
the problem of the complete filling of lakes and fjords which would occur particularly 
when P a is small, but does not solve the problem of the final structure depending 
somewhat on the initial ordering of the perimeter sites. Physically, no recursion rule 
might correspond to the assumption that the formation of additional bonds during 
aging and drying is largely simultaneous and not sequential. 

III. Numerical Results 

Our numerical results for n(X) for the density of states and 7r(A) (related to the 
velocity autocorrelation function of the random walk by a Laplace transform ||, ||) 
are plotted in Fig.2 and Fig. 3 respectively. We see that generally d = 3 data appear 
to be behaving in a clearer way than those of d — 2. We summarize the exponent 
estimates for d s and d w in Tables I and II where the error estimates are mostly from 
the least squares regression and do not take into account any finite size effects or 
other systematic errors that may be present. 
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The grid sizes for the calculations shown in the figures are L = 100 for the square 
lattice and L = 30 for the simple cubic lattice (although we have also checked for 
the cluster size dependence using smaller grids of L = 50 for the square lattice and 
L = 20 for the simple cubic lattice). The spanning and wrapping percolation clusters 
(at p = 0.593 ~ p c for the square and p = 0.312 for the simple cubic lattice) contained 
about 4600 sites for d = 2 and 4000 sites for d = 3 before adding any loops. For 
P Q = 12, the size increases to about 5050 and 4400, for P = 8, to about 5250 and 
4900, and for P = 6, to 6100 and 5700 for d = 2 and d = 3, respectively. For each 
value of P a , averages were taken over 400 independent cluster realizations for the 
square lattice, and for the simple cubic lattice, the number of clusters used varied 
slightly as 481, 480, and 479 for P a = 12, 8, and 6, respectively. 

Thus, both the size of the clusters and the number of independent samples in the 
disorder ensemble are rather modest; however, they are comparable to the previous 
calculations by the same method for the P = oo case || where the results for p c 
were in excellent agreement with other calculations of d s and d w as well as with the 
scaling relation d s = 2df/d w || (but see also |jP2]l). The scope of the calculations 



were CPU time limited mainly because of the need to calculate the eigenvalues and 
eigenvectors very accurately (usually to 6 digits) for the largest 200 or so eigenvalues. 
These calculations took typically 3 hours of CPU time on one Cray Y-MP processor 
for each P Q for the square lattice and 6 hours for each P Q for the simple cubic lattice. 
Performing these calculations on workstations for long periods of time would solve 
the CPU limitation, but then the memory becomes the limiting factor. 

Since the slope in the log-log plot of Fig. 2 must equal d s /2 — 1 (cf. Eq.(fjD), it 
is obvious that the original percolation clusters (P = oo) yield completely different 
values of d s in both d = 2 and 3 even from the case of P Q = 12. Indeed, the estimate 
of d s is about 1.30 for both d = 2 and d = 3 if P Q = oo, but for the loop-added cases, 
the estimates are d s ~ 1.7 for d = 2 and ~ 2.0 for d = 3. Since the case P Q = 12 
contains only about 10 % more sites and visually looks very similar to the original 
percolation cluster (cf. Fig.l (a) and (b)), this is a clear evidence that the structural 
properties are very sensitive to the addition of even a few loops. On the other hand, 
the slopes of the data points for the three finite values of P do not differ very much 
from each other. Thus once the larger loops are added, any further addition of smaller 
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loops do not appear to affect the results very much. 

There does seem to be a tendency for the slope to get steeper when more loops 
are added. This tendency is within the statistical standard deviations for d = 2 but 
definitely outside reasonable statistical errors for d — 3. A study of smaller clusters 
(L = 50 for the square and L = 20 for the simple cubic lattices) also shows different 
behaviors for d = 2 and 3. For the square lattice, this slope gets consistently steeper 
for all P Q = 12, 8, and 6 as L increases; however, the increase in the slope is well 
within the statistical standard deviation in all cases. For the simple cubic lattice, 
on the other hand, the slope becomes less steep for P Q = 12 and steeper for P Q = 6 
as L increases, the changes being slightly outside of the statistical errors. These 
considerations suggest that the different d s for different P Q < oo might be a real effect 
for d — 3 but less likely to be so for d = 2. However, a more complete study, e.g., 
of the finite size effects is needed to definitively answer this question. In any case, 
the differences in the slope among the P Q = 12, 8, and 6 are far smaller than the 
difference between these values of P Q and P Q = oo. 

Compared to the results for the density of states n(A), the results for 7r(A) in 
Fig. 3 show a similar, but much less drastic relative differences in d w between the 
loop-added clusters and the original percolation clusters. That is, the exponent d w 
changes from about 2.9 for the original cluster to about 2.4 for the loop-added case 
for d — 2, a change of only about 17 % compared to the case of d s where the change 
is more than 30 %. This can be understood since the addition of a few loops may 
not change the velocity autocorrelation of a random walker in t steps as much as its 
chances to return to the starting point, or put another way, the structural rigidity 
is much more sensitive to the addition of a few loops than the overall random walk 
displacement on the same network. 

Again, the three-dimensional results appear to be better behaving than the two- 
dimensional ones. However, in both cases, the slopes in the log-log plot of Fig. 3 are 
significantly steeper for P Q = oo than for the loop-added clusters, with a hint of a 
slight trend for the larger slope for P Q = 6 in d = 3. Since this slope must equal 
1 — 2/d w (cf. Eq.(|)), a larger slope implies a larger value of d w . These differences 
among P Q < oo are, however, at least as ill-defined as for the case of n(A). 
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IV. Summary 

In summary, we have presented the analysis of the dynamical properties of the 
random walk confined to loop-enhanced critical percolation clusters. Through the 
vibration-diffusion analogy 0, this work has implications for the (scalar) elastic- 
ity problem on the same fractal network. From the numerical calculations of the 
eigenspectrum of the transition probability matrix, we find what appears to be the 
crossover to a new universality class characterized by a significantly larger spectral 
dimension d s and significantly smaller walk dimension d w as soon as the floppy ends 
of a critical percolation cluster are connected together to form relatively large loops. 
The additional inclusion of successively smaller loops does not seem to further change 
the exponents d s and d w very much. This result supports the observation of Ref . || 
that the deviation between the experimental observation of d s of an aerogel and the 
previous theoretical calculations for the vector elasticity of the critical percolation 
cluster is due to the lack of floppy endings in the aerogel because of the chemical 
reactions during aging and (supercritical) drying. 

The loop-enhanced percolation problem is reminiscent of the problem of the ex- 
ternal surface of the critical percolation cluster. The latter can be defined in several 
ways, e.g., by the percolation hull and by the perimeter sites accessible to a random 
walk from the outside. The accessibility in that case depends on the size of the walker 
|i"3| ; in a somewhat similar way, the floppiness of the loops in our problem depends on 
the minimum loop size P Q . As in our case, the external perimeter problem also leads 
to one new universality class (in addition to the one for the hull) in two dimensions 
fy| . One difference is, however, that our problem leads to a new universality class 
also in three dimensions while this is not so for the external perimeter problem. 

Beyond the qualitative analogy, there might lurk a more quantitative one, since the 
problem of the external surface is that of successively controlling the accessibility for 
a particle diffusing from outside while our problem is that of successively controlling 
the reach of a particle diffusing on the cluster itself. This is, however, beyond the 
scope of this work and must await further research. 
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FIGURES 



Fig. 1: The result of adding successively smaller loops to the critical percolation 
cluster is illustrated with: (a) a typical spanning and wrapping critical percola- 
tion cluster on the square lattice; (b) when all loops of loop parameter P > 12 
are added; (c) when all loops of P > 8 are added; (d) when all loops of P > 6 
are added. The grid size is 100 x 100. 

Fig. 2: Density of eigenvalues n(X) for (a) the square lattice and (b) the simple 
cubic lattice. The symbols x, 0, □, and A correspond to P Q = oo, 12, 8, and 
6, respectively. Lines drawn are linear least squares fits to the corresponding 
data shown here. The data for P a = oo are from Ref. M. 

Fig. 3: The function 7r(A) for (a) the square lattice and (b) the simple cubic lattice. 
The symbols x, 0> n > an d A correspond to P a = oo, 12, 8, and 6, respectively. 
Lines drawn are linear least squares fits to the corresponding data shown here, 
and the data for P Q = oo are from Ref. |J 
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TABLES 



Table 1: Numerical estimates of the exponents d s and d w for loop-enhanced critical 
percolation cluster on the square lattice. The results for P Q = oo are from Ref.8. The 
error estimates are mainly from the least squares regression. 



P d s dy 



oo 1.30 ±0.02 2.86 ±0.04 

12 1.69 ±0.12 2.40 ±0.12 

8 1.66 ±0.10 2.42 ±0.12 

6 1.63 ±0.11 2.39 ±0.15 



Table 2: Numerical estimates of the exponents d s and d w for loop-enhanced critical 
percolation cluster on the simple cubic lattice. The results for P a = oo are from Ref.8. 
The error estimates are mainly from the least squares regression. 



P d s d w 

oo 1.30 ±0.02 3.70 ±0.07 

12 1.99 ±0.02 2.67 ±0.07 

8 1.94 ±0.02 2.76 ±0.08 

6 1.89 ±0.01 2.76 ±0.04 
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